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Abstract 

A recently developed formulation for treating two- and three-nucleon bound states in a three- 
dimensional formulation based on spin-momentum operators is extended to nucleon-nucleon scat- 
tering. Here the nucleon-nucleon t-matrix is represented by six spin-momentum operators accom- 
panied by six scalar functions of momentum vectors. We present the formulation and provide 
numerical examples for the deuteron and nucleon-nucleon scattering observables. A comparison to 
results from a standard partial wave decomposition establishes the reliability of this new formula- 
tion. 

PACS numbers: 21.45.-v, 21.30.-x, 21.45.Bc 



1 



I. INTRODUCTION 



A standard way to obtain scattering observables for nucleon-nucleon (NN) scattering is 
to solve the Schrodinger equation either in momentum or coordinate space by taking advan- 
tage of rotational invariance and introduce a partial wave basis. This is a well established 
procedure and has at low energies (below the pion production threshold) a clear physical 
meaning. At higher energies the number of partial waves needed to obtain converged results 
increases, and approaches based on a direct evaluation of the scattering equation in terms 
of vector variables become more appealing. 

Especially the experience in three- and four-nucleon calculations [H, [2| shows that the 
standard treatment based on a partial wave projected momentum space basis is quite suc- 
cessful at lower energies, but becomes increasingly more tedious with increasing energy, since 
each building block requires extended algebra and intricate numerical realizations. On the 
other hand for a system of three bosons interacting via scalar forces the relative ease with 
which a three-body bound state [3| as well as three-body scattering 0| can be calculated 
in the Faddeev scheme when avoiding an angular momentum decomposition altogether has 
been successfully demonstrated. Thus it is only natural to strive for solving the three nucleon 
(3N) Faddeev equations in a similar fashion. 

Recently we proposed a three-dimensional (3D) formulation of the Faddeev equations 
for 3N bound states [ld\ and 3N scattering ll[ in which the spin-momentum operators are 



evaluated analytically, leaving the Faddeev equations as a finite set of coupled equations 
for scalar functions depending only on vector momenta. One of the basic foundations of 
this formulation rests on the fact that the most general form of the NN interaction can 
only depend on six linearly independent spin-momentum operators, which in turn dictate 
the form of the NN bound and scattering state. Here we extend the formulation of the NN 
bound state given in pj| to NN scattering and provide a numerical realization. 

There have been several approaches of formulating NN scattering without employing a 
partial wave decomposition. A helicity formulation related to the total NN spin was proposed 
in [HI, which was extended to 3N bound state calculations in [6|. The spectator equation 
for relativistic NN scattering has been successfully solved in 0] using a helicity formulation. 
Aside from NN scattering, 3D formulations for the scattering of pions off nucleons (sf and 
protons off light nuclei [9| have recently been successfully carried out. 

In Section [ITl we introduce the formal structure of our approach starting from the most 
general form of the NN potential. We derive the resulting Lippmann-Schwinger equation 
and show how to extract Wolfenstein parameters and NN scattering observables. Numerical 
realizations of our approach that employ a recent chiral next-to-next-leading order (NNLO) 
NN force [iMi as well as the standard one-boson-exchange potential Bonn B [15| are 



presented in Section IIIII The scalar functions, which result from the evaluation of the spin- 
momentum operators and have to be calculated only once are given in Appendices \M and 
[Bl Finally we conclude in Section IIVI The more technical information necessary to perform 
calculations with the chiral potential is given in Appendix O In Appendix [D] the Bonn B 
potential is presented in the form required by our formulation. 

II. THE FORMAL STRUCTURE 

We start by projecting the NN potential on the NN isospin states | tnit), with t = 0,mt = 
being the singlet and t = l,mt = — 1,0,1 the triplet. We assume that isospin is conserved, 
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but allow for charge independence and charge symmetry breaking, and thus for a dependence 
on rrit, 

{t'm[\V\tmt) = Su'Sm,m'y''^' (2.1) 
Furthermore, the most general rotational, parity and time reversal invariant form of the 



off-shell NN force can be expanded into six scalar spin- momentum operators [17|, which we 
choose as 

Wi{(Ti,(T2,p',p) = 1 

W2i(Ti,(T2,p',p) = (Ti - (T2 

WsicTl, (T2, p', p) = « (CTI + 0-2) • (p X p') 

Wi{(Ti, a2, p', p) = CTi ■ (p X p') 0-2 • (p X p') 

W^{cri, (T2, p', p) = CTi ■ (p' + p) 0-2 ■ (p' + p) 

Wq{(Ti, era, p', p) = cr^ • (p' - p) 0-2 ■ (p' - p) (2.2) 

Each of these operators is multiplied with scalar functions which depend only on the mo- 
menta p and p', leading to the most general expansion for any NN potential 

6 

ytm. ^ J2v^m.{p',p) W,{cT,, CT2 , p' , p) (2.3) 

i=i 

The property of Eq. (12. ip carries over to the NN t-operator, which fulfills the Lippmann- 
Schwinger (LS) equation 

with Go{z) = {z — Hq)~^ being the free resolvent. The t-matrix element has an expansion 
analogous to the potential, 

6 

^tm. ^ J2t'p{p',p) w,{a^,a2,p',p) (2.5) 

i=i 

Inserting Eqs. (12. 3p and (12. 5 p into the LS equation (12. 4p . operating with Wk{(Ti, cr2,p' ,p) 
from the left and performing the trace in the NN spin space leads to 

E^'^^(p''P)^r(p''P) = E^^'^(p''P)^r(p''P) 

j 3 

+ / d?p"Y,v'p{p\p")G,{p") t;.r(p",p) S.,y(p',p",p). (2.6) 

33' 

The scalar coefficients Ai^j and Sfcjj' are defined as 

Akj{p',p) = Tr(^zi;fc(cri,cr2,p',p) WJj(cri,cr2,p',p)j (2.7) 

Bkjj'{p\p'\p) = TT(wk{(Ti,(T2,p',p) Wj{(Ti,(T2,p',p") Wj>{(Ti,(T2,p",p)^ (2.8) 

Here all spin dependencies are analytically evaluated, and the coefficients only depend on the 
vectors p, p', and p". The explicit expressions for the coefficients are given in Appendix El 
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Thus we end up with a set of six coupled equations for the scalar functions i*™'*(p', p), 
which depend for fixed |p| on two other variables, |p'| and the cosine of the relative angle 
between the vectors p' and p, given by p' ■ p. 

Since Eqs. fl2.3p and f l2.5p are completely general, any arbitrary NN force can be cast into 
this form and serve as input. Finalizing the formulation, we only need to antisymmetrized 
in the initial state by applying (1 — Pi2)\p)\mim2)\t'mt) , and consider the on-shell t-matrix 
element for given trrit: 



on— shell 



HI 

M*T , = (27r)2t*T , (p',p) 

= -^(27r)2(Km'2|[t*'"*(p',p) + (-)*t*™*(p',-p)P^2]l"^i"^2)). (2.9) 

Here PI2 interchanges the spin magnetic quantum numbers for the initial particles, m rep- 
resents the nucleon mass. 

For the on-shell condition, characterized by |p'| = |p| the vectors p — p' and p + p' 
are orthogonal. Under this condition, the operator cri ■ 02 can be represented as a linear 
combination of the operators Wj^ j = A — 6 [isj, i.e. 

CTl ■ 0-2 = ^— ^ CTi • (p X p') 0-2 • (p X p') 

(P X P )^ 

+ (p^p/)2 '^l ■ (P + P') ^^2 ■ (P + P') 

+ (p_^p.)2 '^l ■ (P - P') • (P - P') (2.10) 

We can use the relation of Eq. fl2.10p for internal consistency checks of the calculations. 
However, in order to keep the most general off-shell structure of Eq. ( 12. Sp . we need to keep 
all six terms. We will come back to the numerical implications of this fact below. 
From Eq. f l2.9p we read off that the scattering matrix is given by 
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m 



M'r , = --(27r)2 y 

+ (-)*4"*(P'' -P) {m[m2\wj{ai,a2,p',-p)\m2mi) 



t^Pip', P) {m[m2\wj{(Ti, 0-2, p', p)|mim2) 



(2.11) 



On the other hand the standard form of the on-shell t-matrix for given quantum numbers 
trrit [l8| reads in the Wolfenstein representation 

M*T , = a*"* {m\m'\mim2) 

- i-, 7| {m[m'2\w3{(Ti,(T2,p',p)\mim2) 

IP X p I 

+ I -^{m[m2\w4{cri,(T2,p',p)\mim2) 

IP X p I 

+ 7 — \ — -IT {m[m2\w5{(Ti,(T2,p',p)\mim2) 
(P + pr 

+ J {m[m'2\wQ{cri,(T2,p',p)\mim2) (2.12) 
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Due to the action of Eq. (12. 9p . which interchanges mi with 1122, the two parts 

of Eq. (12. lip yield different results. Again, standard relations [l8|, ll9| must be applied to 
extract the Wolfenstein parameters: 



- Tr (M) 
-i- Tr I M 



m 



tint 



- Tr ( M 



|P X P'l 
w^4(cri,cr2,p',p) 



/|2 



{9 + h)'^^ 
{g - hf^^ 



- Tr ( M 

- Tr ( M 



|P X P 
w^5(cri,cr2,p',p) 



'^2 



(P + P 
w^6(cri,cr2,p',p) 



IP-P 



'^2 



(2.13) 



It is straightforward to work out those relations starting from Eq. (12. lip . In order to simplify 



the notation we write tj = tf'ip', p), tj = t*'"'(p', 



-p), and X = p' ■ p and obtain 



a 



tmt 



tl + (-) 

ip" 



1_3_1 



[h - (-)%) 



m 



tmt 



\h - \i2 + ^-p\l - x^)h 



tmt 
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j^tm 



-p'{l-x%-p'{l + x% 



t2 +p\l + X)h + p'(l - X)U + (-) 



\h - \i2 - \p\l - x'Yu 



p^il + x)h - - x)U + (-)* - - x% + p2(l + x)t 



(2.14) 



It remains to consider the particle representation. For the proton-proton or neutron-neutron 
system the isospin is t = 1. Thus the above given Wolfenstein parameters are already the 
physical ones and enter the calculation of observables. In the case of the neutron-proton 
system both isospins contribute and the physical amplitudes are given by \{aP^ + a^°), 
i(cOO + c^O),etc. 

Once the Wolfenstein parameters are knowri, all NN observables can readily be calculated 
taking well defined bilinear products thereof . For example, the spin averaged differential 
cross section /q is given as |Tr MM^ . 

For completeness, we also give the derivation of the deuteron which carries isospin t = 
and total spin s = 1. We employ the operator form from Ref. [5|, 



(Pl^m.) 



01 (P) + ( O-l • P 0-2 • P - ) 02(P) 



2 

E 

k=l 



(f>k{p) hk{(Ti,a2,p)\lmd), 



(2.15) 



where \lmd) describes the state in which the two spin-| states are coupled to the total spin-1 
and the magnetic quantum number m^i. The definition of the operators hk can be easily read 
off the first line of Eq. (I2.15p . The two scalar functions 0i(]3) and 02 (p) are related in a 
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simple way to the standard s- and rf-wave components of the deuteron wave function, i/jqIp) 
and V'2(p) by [5| 



Mp) = ^'^2b)- 

Next we use the Schrodinger equation in integral form projected on isospin states, 



(2.16) 



(2.17) 



Inserting the explicit expression of Eq. (12.151) we obtain 

lip) + f CTl ■ P ■ P - ^PM 02(P) 



X 



1 f 

_p2 / dVXl^r(P'P) ^i('^l''^2,P,p') 

Mp') + f ■ p' 0-2 ■ p' - ^p'^j Mp) l^ma), 



(2.18) 



where is the deuteron binding energy. We remove the spin dependence by projecting 
from the left with {lmd\bk{(Ti, (T2, p) and summing over nid- This leads to 

1 2 

{lmd\bk{(Ti,a2,p) ^(j)k'ip)bk'icri,cr2,p)\lmd) = 



md=-l 



Ed 



k'=l 
6 



^ J a?VX^^^f(P,p')«'j(o-i,o-2,P,p') 

m md=-l j=l 



^ (pk'ip') 6fc"(cri,cr2,p')|lmd) . 



(2.19) 



fc"=i 



Defining the scalar functions 



^kk'iP)= Y {lmd\bk{cTi,CT2,p)bk'{(Ti,a2,p)\lmd) 

md=~l 



(2.20) 



and 



Btjk"iP^P')= X] {^md\bkicri,cT2,p)wj{(Ti,a2,p,p')bk'i(Ti,(T2,p')\lmd), (2.21) 

md=-l 

we obtain for Eq. fl2.19p 

E ^iMMiP) = / rfV E ^r(P' p') E ^'.^"(P' p')M'iP') . (2.22) 



fc'=l 



k"=l 



Note that A^^, and B^- are both independent of the interaction. Therefore, these coeffi- 
cients can be prepared beforehand for all calculations of the deuteron bound state, which 
consists of two coupled equations for the functions and (j)2{j>)- The summation over 

md guarantees the scalar nature of the functions and Bfji^„{p, p'), which are given in 

Appendix [HI The azimuthal angle can be trivially integrated out, leading to the final form 
of the deuteron equation 



dp'p"Mp') dxY,vf{p.p\x)B%AP.P\x), (2.23) 

m k"=l 1=1 



k'=l 



where x = p' ■ p. 

III. NUMERICAL REALIZATION 
A. The deuteron 

For a numerical treatment of Eq. (12.231) , it is convenient to first define 



/I 6 



(3.1) 



and then assume that the integral over p' will be carried out with some choice of Gaussian 
points and weights {pj,gj) with j = 1, 2, . . . , A^. This leads to 



[9oP]Zkk'{Pi.P,) + 5,j^Ai^,{p,^Uy{p^) = EdY^^tk'iP^)^k'{p^)■ (3.2) 

k'=l j=l ^ ^ k'=l 

Eq. (13.21) can be written as a so-called generalized eigenvalue problem 

= E,Y^, (3.3) 



or 



2N 2N 



YRii'^i' = EdY^u<i', (3.4) 



i'=i i'=i 



where 



/ = ^ + (A; - i)Ar, 
6' = <Pk'iPj), l'=j + {k'-l)N, 



p ■ 

Rw = 9jP]Zkk'{Pi,Pj) + 5ij^^A'lj^,{pi) 
1 



2m7r 



yw = ^^l^Kk'iPi)- (3.5) 



TABLE I: The parameters of the chiral potential of Ref. 13| in order NNLO. The LEC's are given 
for the cutoff combination A= 600 MeV and A= 700 MeV. The pion decay constant and masses 
are given in MeV. The constants Cj are given in GeV~^, Cs and Ct in GeV~^ and the other Ci in 
GeV-4 



9A 
1.29 


92.4 


134.977 


139.570 


m 
938.919 


Cl 

-0.81 


C3 
-3.40 


C4 

3.40 




Cs 
-112.932 


2.60161 


Ci 
385.633 


1343.49 


C3 
-121.543 


C4 
-614.322 


C5 
1269.04 


Ce 
-26.4880 


Cr 
-1385.12 



TABLE IL Meson parameters for the Bonn B potential 15|]. The a parameters shown in the table 
are for NN total isospin 0. For NN total isospin 1 they should be replaced by mg- = 550 MeV, 
^.9437, A„ = 1.9 GeV and n = l. 



9a, 

An 



meson 


rtia [MeV] 


ai 

47r 




A. [GeV] 


n 


vr 


138.03 


14.4 




1.7 


1 


V 


548.8 


3 




1.5 


1 


6 


983 


2.488 




2 


1 


a 


720 


18.3773 




2 


1 


P 


769 


0.9 


6.1 


1.85 


2 




782.6 


24.5 





1.85 


2 



Since A'^^ 7^ 0, = ^21 ~ ^^"^ ^22 7^ 0' matrix Y is diagonal and can be easily 
inverted, we encounter an eigenvalue problem 

[Y-'R) C = E,^, (3.6) 

which is of the same type and dimension as is being solved for the deuteron wave function in 
a standard partial wave representation, where one calculates the s- and ci-wave components, 
ipo{p) and ip2{p)- The connection between the two solutions, {(f>i{p), 02 (p)) and ("00 (p)) "02 (p)), 
given by Eqs. (12.161) provides a direct check of the numerical accuracy. 

As a first example we use a chiral NNLO potential which for the convenience of 
the reader is briefly described in Appendix C. For the specific calculation performed here 
we take the neutron-proton version of this potential and employ the parameters listed in 
Table [11 

We consistently use these potential parameters in the 3D and the PW calculations. In 
the first case we solve Eq. (13.61) for 0i(p) and 02 (p) and then use Eqs. (I2.16P to obtain ipoip) 
and ip2{p)- In the second case we employ the standard partial wave representation of the 
potential and solve the Schrodinger equation directly for ilJoip) and "02 (p)- Both methods 
give the same value for the deuteron binding energy, namely £"^=-2. 19993 MeV and s-state 
probability Ps=95.291 %. The wave functions are identical as can be seen in Fig. [T] 

As second NN force we choose the Bonn B potential [l5|, which has a more intricate 
structure due to the different meson-exchanges and the Dirac spinors. The operator form 
of this potential, corresponding to the basis of Eq. (12.21) is derived in Appendix [D] and the 
parameters are given in Table HTl In this case the nucleon mass is set to m= 939.039 MeV. 
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Again we have an excellent agreement between the 3D and the partial wave based calculation 
for the deuteron binding energy, £^^=-2.2242 MeV, the s-state probability {Ps= 95.014 %) 
and the wave functions, which are displayed in Fig. |2l 

In summary, we confirm that the 3D approach gives numerically stable results, which are 
in perfect agreement with the calculations based on standard partial wave methods. 




0.5 1 1.5 2 1 2 3 4 

p[fm-^] p[fm-^] 



FIG. 1: The s-wave (left) and d-wave (right) component of the deuteron wave function as a function 
of the relative momentum p for the chiral NNLO potential specified in the text. Crosses represent 
results obtained with the operator approach and solid lines are from the standard partial wave 
decomposition. 



B. NN scattering observables 

The inhomogeneous LS equation fl2.6p for the six components t*™' can be solved for a 
fixed value of p. For the vectors p and p' we choose the explicit representation 

p = (0,0,1) 



1 -x'^0,2;') 



p" = {VI- x"'^cosip",Vl- x"^smip",x") (3.7) 

so that the scalar products become 

p' ■ p = x' 
p" ■ p = x" 

p' ■ p" = x'x" + Vl - x'Wl - x"^cos if" = y. (3.8) 



Let us now calculate the integral term on the right-hand-side of Eq. f l2.6p for a positive 
energy of the NN system, E'c.m 



m ' 



P 



Skip',p,x') ^ I dfp"'- jiT-- fkip";p',P,x'), (3.9) 
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0.5 1 1.5 2 1 2 3 4 

P [fm"'] p [fm-^] 



FIG. 2: The same as in Fig. [T]but for the Bonn B potential [15]. 



where 



27r 



J J B,,,ip',p",p,x',x",ip") v,{p',p",y) tAp",p,x") . (3.10) 
j,i'=i „i 

Here the index trrit for the t-matrix element is omitted for simphcity. For the momentum 
integration in Eq. (13. 9p an upper bound p is introduced, since the contributions to the 
integral for larger momenta are insignificant, the potential and the t-matrix are essentially 
zero. Then the integral of Eq. (13. 9 p can be treated in a standard fashion and one obtains 

p 

Sk{p,p,x) = / dp ^ -2 + -Pofkipo) In^ m . 3.11 

J Po-P 2 \ P-Po J 



It is tempting to solve Eq. (12. 6p by iteration and then sum the resulting Neumann series 
with a Fade scheme. The determinant of the 6x6 matrix A{p' ,p,x'), which appears on 
both sides of (12. 6p . can be easily calculated with the result 

det{A) = -65536 j9>'* (p^ -p'^)^ (l - x'^Y . (3.12) 

In particular, this determinant is zero for p' = p and x' = ±1. However, by a careful choice 
of the p, p' and x' points, it is possible to work with non-zero values of det(74), so that the 
matrix A can be inverted. In this case Eq. (12.61) can be written as 

t{p',p,x') = v{p',p,x') + A-^{p',p,x') S{p',p,x'), (3.13) 

where t{p',p,x'), v{p',p,x') and S{p',p,x') denote now six-dimensional vectors with compo- 
nents tj, Vj and Sj. Note that S{p',p,x') contains the unknown vector t{p',p,x'). We arrive 
at the following iteration scheme: 

t''^\p' ,p, x') = v{p',p,x') 

&\p',p,x') = v{p',p,x') + A-\p',p,x')S^''-^\p',p,x'), forri>l, (3.14) 
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where S^'^~^\p',p,x') is calculated using the vector t{p',p,x') from the previous iteration, 
i.e. t^"'~^^{p',p,x'). However, our experience with this iteration scheme is discouraging. 
Numerically det(v4) can be very close to zero, and in such cases the rank of matrix A can 
vary from 2 to 5. As a consequence, it is very difficult to maintain numerical stability for 
this iterative method. Another drawback of using the inverse of A is that it is impossible to 
obtain the on-shell matrix element t{po,po, x') directly. One would have to rely on numerical 
interpolations for calculating on-shell matrix elements. 

For this reason we decided to solve Eq. (12. 6p directly as a system of inhomogeneous 
coupled algebraic equations. To this aim we first perform a discretization with respect to 
the different variables in the problem. As typical grid sizes we take Ux = 36 Gaussian points 
for the x" integration, and use the same grid for the x' points. Furthermore, we use Up = 36 
Gaussian points for the p' and p" grids, which are defined the interval {0,p = 40fm~^). 
These points are distributed in such a way that po is avoided and the same number of points 
is put symmetrically into two narrow intervals on each side of po 0]. Such a choice proved 
advantageous in the treatment of the ^Sq channel for the PWD calculations and is kept here. 
In addition, po is added to the set of p' points. Finally, we choose n^// = 60 Gaussian points 
for the ip" integration. Thus, we arrive at a system of 6 x (n^ + 1) x linear equations of 
the form 

= b, (3.15) 

where the vector ^ represents all unknown values of tj{p',p, x') for fixed p. If we choose from 
the very beginning p = Po, then the solution of Eq. (I3.15P contains the on-shell t-matrix in 
the operator form, namely tj{po,po, x'). 

It is clear that for the on-shell t-matrix the solution cannot be unique, since the six oper- 
ators become linearly dependent on each other (see Eq. (12.101) ). In principle, one therefore 
expects that Eq. (13.151) is non- invert ible and that tools like a singular value decomposition 
are required for the solution. However, we found that this is not required since the stan- 
dard LU decomposition of Numerical Recipes worked safely for both interactions, all 
the considered laboratory energies and different choices of the mesh points. Interestingly, 
the actual solution for the on-shell t-matrix is not unique as expected and depends even on 
the optimization level of the compiler. However, the observables turn out to be stable and 
unique. 

Of course, setting p = pq is not necessary. For p ^ po the system of equations (I3.15P has 
a unique and smooth solution and afterwards the interpolation to the on shell case can be 
safely performed. 

The path to NN observables is straightforward. From Eq. (12.111) we evaluate first the 
scattering matrix M for all possible spin projections m'^, m'g, mi, and m2, noting that 
on-shell 

tr(p''P)=^r(po,Po,x') (3.16) 

and 

t'p{v'.-p)=t'p{p,,Po-x'). (3.17) 

Since we use a set of x' points which is symmetric with respect to z' = 0, no interpolation 
is required and M is easily obtained. Before we can make use of Eq. (I2.13p . we calculate 
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matrix elements of the modified operators wj appearing in fl2.13p . in the same representation 
as for the matrix M: 



m'^ 777-2 



777'^ 7772 



m\m'2 



m[m'2 







p) 




|P X P'l 




104 


(cri,cr2,p' 


P) 




P ^ P' P 




W5 


(cri,cr2,p' 


P) 




(P + P')' 




Wei 


Ci, cr2, p', 


P) 


(P - P')' 



777 i7772 
777 i 7772 
777 i 7772 
777 i 7772 



(3.18) 



For this calculation symbolic software like Mathematica® 2l| proves very useful. In the 
next step, the Wolfenstein parameters are calculated as sums over 777'^, 7772, 7771 and 7772. For 
example 



777x7772 



«'3lcri,cr2,P ,Pj 



|P X p'l 



m\m'2 



(3.19) 



Finally, the NN observables result from the Wolfenstein parameters as simple bilinear ex- 
pressions 
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In Figs. ISHHlwe compare a selected set of observables calculated with the new 3D method 
to results obtained by using a standard partial wave decomposition, employing the same 
potentials we used for the deuteron calculations. For the chiral potential we chose two 
laboratory kinetic energies 13 and 150 MeV, whereas for the Bonn B potential the higher 
energy is chosen to be 300 MeV. We made sure that in all cases a sufficient number of partial 
waves is included to obtain converged results in the standard PWD approach. For all the 
energies considered our converged PWD results agree perfectly with predictions obtained 
from the new 3D approach. 

In Figs. [THE] we demonstrate the convergence with respect to different maximum total 
angular momenta jmax towards the results calculated using our new 3D method for the 
differential cross section and the asymmetry A. Here we employ the Bonn B potential and 
show the calculations for the neutron-neutron and neutron-proton cases separately. As one 
can see, quite a sizeable number of partial waves is required for a converged calculation 
at 300 MeV. Finally, in Fig. |9] we display the Wolfenstein amplitudes for neutron-proton 
scattering at 300 MeV laboratory kinetic energy. Again we compare partial wave based 
calculations for different maximum total angular momenta jmax to the 3D calculation. We 
observe that the maximum number of partial waves needed for obtaining a converged result 
is quite different for the different amplitudes. 



IV. SUMMARY AND CONCLUSIONS 

Two nucleon scattering at intermediate energies of a few hundred MeV requires quite 
a few angular momentum states in order to achieve convergence of e.g. scattering observ- 
ables. We formulated and numerically illustrated an approach to treating the NN system 
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FIG. 3: Selected observables for the neutron-neutron (left panel) and neutron-proton (right panel) 
system at the projectile laboratory kinetic energy 13 MeV as a function of the center of mass angle 
6 for the chiral NNLO potential 1^. Crosses represent results obtained with the operator approach 
and solid lines represent fully converged results from the standard PWD. For the definition of the 
R, A and D observables see e.g. 18l |. 
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FIG. 5: The same as in Fig. [3] for the projectile laboratory kinetic energy being 150 MeV. 
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FIG. 6: The same as in Fig. Hlfor the projectile laboratory kinetic energy being 300 MeV. 
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FIG. 7: The convergence of the PWD results for the differential cross section and the depolariza- 
tion coefficient A for neutron-neutron scattering based on different numbers of partial waves 
determined by the maximal total angular momentum jmax of the NN system (lines) with respect 
to the result of the three-dimensional calculation (crosses) for projectile laboratory kinetic energy 
300 MeV and the Bonn B potential [15|]. 

working directly with momentum vectors and using spin-momentum operators multiplied 
by scalar functions, which only depend on the momentum vectors. This approach is quite 
natural, since any general NN force being invariant under time- reversal, parity and Galilei 
(or Lorentz) transformations can only depend on six linear independent spin-momentum op- 
erators. The representation of the NN potential using spin-momentum operators leads to a 
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FIG. 8: The same as in Fig. [7] for neutron-proton scattering. 



system of six coupled equations of scalar functions (depending on momentum vectors) for the 
NN t-matrix, once the spin-momentum operators are analytically calculated by performing 
suitable trace operations. 

We calculated deuteron properties and NN scattering observables using two different NN 
potentials, one derived from chiral effective field theory and one from meson exchange. For 
all cases we found perfect agreement between the calculations based on our new method and 
conventional calculations using a partial wave basis. 

This work is intended to serve as starting point towards treating three-nucleon systems 
without partial waves. The theoretical formulation has already been given for the 3N bound 
state (including 3N forces) in 10 1 and for 3N scattering in [ll|. For a much simpler case when 
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FIG. 9: The Wolfenstein parameters for neutron-proton scattering for projectile laboratory kinetic 
energy 300 MeV calculated with the Bonn B potential 0. Results of the 3D calculatfon are given 
by the crosses. The convergence of the PWD results for increasing values of maximum angular 
momentum jmax is shown by the different curves labeled in the figure. The left panels show the 
real parts of the amplitudes, whereas the imaginary parts are displayed in the right panels. 
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spin- and isospin-degrees are neglected the feasibility of three-body scattering calculation in 
the Ge V regime has already been demonstrated 0] , even including Poincare symmetry 22 



Since our approach leads to coupled equations of scalar functions of momentum vectors, the 
generalization to include spin-degrees of freedom appears feasible. In the 3N system not only 
the number of partial waves increases rapidly, but also 3N forces appear as new dynamical 
input. In particular in the chiral approach the number of 3NF contributions proliferates 
with order of expansion of the theory. In higher orders many complicated terms contribute 
to the 3N force [l3|. In this case a traditional partial wave decomposition of the 3NF poses a 
serious problem which has a chance to be alleviated by a direct three-dimensional treatment. 
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Appendix A: Coefficients for NN scattering 

In this appendix we present the expressions A and B given in Eqs. f l2.7p and fl2.8p for 
NN scattering. The coefficients Aij{p',p) can be obtained in terms of the following four 
functions FA 

FA,{p',p) = 4(pxp')2 (Al) 

FA(p',p) = 4(p' + p)^ (A2) 

FA,{p',p) = 4(p'-p)2 (A3) 

FA(P',P) = ^ip"-py (A4) 
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The non-zero coefficients Aij{p',p) are: 





= 4 




^22(P',P) 


= 12 




^24(P', P) 


= A42(P',P) = 


i^^3(p',p) 


^25 IP , P) 


= A52(p',p) = 


FA4(p',p) 


^26 (P', P) 


= A62(p',p) = 


i^^(p',p) 


^33(P', P) 


= -2FA3(p',p) 


^44(P', P) 


= 4^24(P',P) 




A5(P',P) 


= |^^5(P',P) 




A6(P',P) 


= A65(p',p) = 


i^A(p',p) 


^66(P',P) 


= |^26(P',P) 





All other Aij{p',p) 



0. 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 

(All) 

(A12) 
(A13) 
(A14) 



The non-vanishing coefficients Bikj{p', p", p) can be expressed by means of the following 
25 functions FB: 



FBsaip', P",P) 
FS3ft(p',p",p) 
F53e(p',p",p) 



4(p X p")' 
4(p" X p')^ 
4(p X p')' 



(A15) 
(A16) 
(A17) 



F54„(p',p",p) 
FB,,{p',p",p) 

FS4e(p',p",p) 



4(p" + P)^ 

4(p' + P")' 
4(p' + pf 



(A18) 
(A19) 
(A20) 



fbUp', p", p) 

FB56(P',P",P) 
FB5C(P',P",P) 



4(p" - P)^ 
4(P' - P")' 
4(p' - P)' 



(A21) 
(A22) 
(A23) 



FB6„(p',p",p) = -8(p" X p') . (p X p") 

FS66(p',p",p) = -8(p X p') ■ (p" X p') 

FS6e(p',p",p) = -8(p X p') • (p X p") 



(A24) 
(A25) 
(A26) 



m(p',p",p) = 4{(pxp')-p'r 



(A27) 



FB8a(p',p",p) = 2{(p' + p")-(p" + p)} 
FB8,(p',p",p) = 2{(p' + p).(p' + p")} 
FS8c(p',p",p) = 2{(p' + p)-(p" + p)} 



(A28) 
(A29) 
(A30) 
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FBUp', p", p) = 2{(p' - p") ■ (p" - p)} (A31) 

FB,,ip', p", p) = 2{(p' - p) ■ (p' - p")} (A32) 

FB9e(p', P", P) = 2{(p' - p) • (p" - p)} (A33) 

FB.oaip'., P", P) = 2{(p' + p") . (p" - p)} (A34) 

FSio.(p', P", P) = 2{(p' + p) . (p' - p")} (A35) 

FSioc(p', p", P) = 2{(p' - p) ■ (p" + p)} (A36) 

FBUp', p", p) = 2{(p' - p") • (p" + p)} (A37) 

FBnbip', P", P) = 2{(p' - p) • (p' + p")} (A38) 

FBndp': P", P) = 2{(p' + p) • (p" - p)} (A39) 

The non-zero Sjfej(p', p", p): 

BMp', p", p) = B2u{p', p", p) = 522i(p', p", p) = 12 (A40) 

5i24(p', p", P) = 5214(P', p", p) = FB,,{p', p", p) (A41) 

5i25(p', p", P) = B2u(p', p", P) = FB^aip', p", p) (A42) 

5l26(p', P", P) = 52i6(p', P", P) = FB,a{p', P", P) (A43) 

5133(P', P", P) = 5233(P', P", P) = FB(,a{p', p" , p) (A44) 

5i42(p', p", p) = B24i(p', p", p) = FBs,(p', p", p) (A45) 

5i44(p',p",p) = ^FB,a{p',p",pf (A46) 

5l45(p', P", P) = Si46(p', P", P) = Si54(p', P", P) = Bl64(p', p", p) 
= 54i5(p', P", P) = S4i6(p', p", p) = B514(P', P", P) 
= S6i4(p', P", P) = S45l(p', p", p) = B46l(p', P", P) 

= S54l(p', p", p) = S64l(p', p", p) = FB,{p', p", p) (A47) 

5i52(p', p", P) = 525i(p', P", P) = FB^,{p', p", p) (A48) 

5i55(p',p",p) = FBsa{p,p",pf (A49) 

5i56(p',p",p) = FB,oa{p',p",P? (A50) 

i?i62(p',p",p) = i?26i(p',p",p) = FB,,{p',p",p) (A51) 

5i65(p',p",p) = FBna{p,p",py (A52) 

5i66(p',p",p) = FBga{p',p",pf (A53) 

Bsisip', p", p) = S323(p', p", P) = FBe,{p', p", p) (A54) 

5412(P', p", P) = 542i(p', p", p) = FBscip', p", p) (A55) 

5414(P',P",P) = -^FB6e(p',p",p)' (A56) 

lb 

55i2(p', p", P) = 552i(p', p", P) = FB,,{p', p", p) (A57) 

5515(P',P",P) = FS8e(p',p",p)' (A58) 

i?5i6(p',p",p) = FSnc(p',p",p)' (A59) 

5612(P', P", P) = B62l(p', P", P) = FB5c(p', P", P) (A60) 
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^615(P ,P ,Pj 


= FSioc(p,p ,p) 


(A61) 




— ^ ^9c\y ) P ) Pj 




R„oi ( r\' rt" rt\ 
-D33ll,P ) P ) PJ 


— ^332\y ,P ) PJ — -f^ -D66l,P ) P ) PJ 




544l(p',P",P) 


= ^FB66(P',P",P)' 


(A64) 


-S551(P',P",P) 


= FS8fe(p',p",p)' 


(A65) 


R (■r^.' r^" 
-D56IIP , P , Pj 


— -f'-DlOfelP >P )Pj 




^65l(p',P",P) 


= FSn6(p',p",p)^ 


(A67) 


^66UP ,P ,Pj 


1—1 -I— » / / // \ 9 

= FS96(p',p",p)2 


(A68) 


Bill Co' d" d"! 


= 4 


(A69) 


5244(P',P",P) 


= -/'FB,(p',p",p) 


(A70) 


5245(P',P",P) 


= -^F536(p',p",p)ma(p',p",p) 






+m(p',p",p) 


(ATI) 


5246(P',P",P) 


= -iFS35(p',p",p)FS5„(p',p",p) 






+FS7(p',p",p) 


(A72) 


i^254(p',p",p) 


= -^FS3a(p',p",p)m;,(p',p",p) 






+FS7(p',p",p) 


(A73) 


5255(P',P",P) 


= -^F545(p',p",p)F54a(p',p",p) 
4 






+FS8„(p',p",p)^ 


(A74) 


5256(P',P",P) 


= -iFi?46(p',p",p)FB5a(p',p",p) 






■1—1 1— » / I II \ 

+FSiOa(p',p",p)' 


(A75) 


5264(P',P",P) 


= -^FS3„(p',p",p)F556(p',p",p) 






+FS7(p',p",p) 


(A76) 


5265(P',P",P) 


= -iFB56(p',p",p)Fi^4a(p',p",p) 






+FSna(p',p",p)' 


(A77) 










+FB,,{p',p",pf 


(A78) 
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5424(p',p",p) = Vm(p',p",p) (A79) 

5425(P',P",P) = -^FS3e(p',p",p)FS4„(p',p",p) 

+Fi?7(p',p",p) (A80) 

5426(P',P",P) = -^Fi?3c(p',p",p)Fi^5a(p',p",p) 

+F57(p',p",p) (A81) 

5524(P',P",P) = -^FS3„(p',p",p)F54e(p',p",p) 

+Fi?7(p',p",p) (A82) 

5525(P',P",P) = -^Fi?4a(p',p",p)mc(p',p",p) 

+FS8c(p',p",p)' (A83) 



5526(P',P",P) - -^FS5„(p',p",p)mc(p',p",p) 

+FSnc(p',p",p)' (A84) 

5624(P',P",P) = -^FB3„(p',p",p)FB5c(p',p",p) 

+FS7(p',p",p) (A85) 

5625(P',P",P) = -^FS4„(p',p",p)F55e(p',p",p) 

+FSioc(p',p",p)' (A86) 

5626(P',P",P) = -^FS5a(p',p",p)mc(p',p",p) 

+F59,(p',p",p)2 (A87) 
5442(p',p",p) = -p''m(p',p",p) (A88) 

5452(P',P",P) = -^FS3e(p',p",p)FS46(p',p",p) 



4 

+m(p',p",p) (A89) 



5462(P',P",P) = -^FS3e(p',p",p)F555(p',p",p) 



+F57(p',p",p) (A90) 

5542(P',P",P) = -^FS3b(p',p",p)FS4c(p',p",p) 

+m(p',p",p) (A91) 
p',p",p) 

+FS8;,(p',p",p)2 (A92) 

3',p",p)i 

+FSio6(p',p",p)' (A93) 



5552(P',P",P) = -^me(p',p",p)F54,(p',p",p) 



5562(P',P",P) = -^me(p',p",p)FS55(p',p",p) 
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5642(P',P",P) = -^FS3b(p',p",p)F55e(p',p",p) 

+m(p',p",p) (A94) 

i?652(p',p",p) = -^FS5c(p',p",p)Fi346(p',p",p) 

+FBnb{p',p",pf (A95) 

5662(P',P",P) = -^FS5e(p',p",p)FS55(p',p",p) 

+FS9(.(P',P",P)' (A96) 

i?224(p',p",p) = -2FS3a(p',p",p) (A97) 

5225(p',p",p) = -2FS4„(p',p",p) (A98) 

S226(P,P",P) = -2FB5a(p,p",p) (A99) 

i?242(p',p",p) = -2FS3fe(p',p",p) (AlOO) 

i?252(p',p",p) = -2F54,(p',p",p) (AlOl) 

5262(P',P",P) = -2Fi?5b(p',p",p) (A102) 

5422(p',p",p) = -2Fi?3c(p',p",p) (A103) 

5522(P',P",P) = -2me(p',p",p) (A104) 

5622(p',p",p) = -2F55e(p',p",p) (A105) 

5222(p',p",p) = -24 (A106) 

i?344(p',p",p) = \FBea{p',p",p)FBr{p',p'\p) (A107) 

5345(p', P", P) = 5435(p', P", P) = -2{(p" + p) ■ p'}FB,{p', p", p) (A108) 

i^346(p', P", P) = 5436(p', P", P) = -2{(p" - p) • p'jFBrip', p", p) (A109) 

5355(P',P",P) = ^{2FS3,(p',p",p)-FS6,(p',p",p) 

-FB,,{p', p", p)} FBsaip', p", p) (Alio) 

i^356(p',p",p) = -^{2FB,,{p',p",p) + FB,,{p', p",p) 

-FBedp', p", p)} Fi3ioa(p', P", P) (Alll) 

5365(P',P",P) = ^{2FS3c(p',p",p)-FS66(p',p",p) 

+FS6e(p', p", p)}FBna{p', p", P) (A112) 

i?366(p',p",p) = -^{2FS3c(p',p",p)+Fi?6f.(p',p",p) 

+FBec{p', p", p)} FB,a{p', p", p) (A113) 

5434(P',P",P) = ^F56e(p',p",p)m(p',p",p) (A114) 
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5534(p', P", p) = S543(p', p", P) = 2{(p' + p) • p"}FBr{p', p", p) (A115) 
5634(p', P", P) = -5643(p', p", P) = 2{(p' - p) . p"}FB,{p', p", p) (A116) 

5535(P',P",P) = -^{2FS3;,(p',p",p) + FS6„(p',p",p) 

-FBebip', p", p)} FBseip', p", p) (A117) 

i?536(p', P", P) = {2FS36(P , P", P) + Fi?6a(p', p", p) 

+FBe,{p', p", p)} FBucip', P", P) (A118) 

5635(p', P", P) = {'2FB,,{p', p", p) - FB,,{p', p", p) 

-FS6b(p', p", p)} FSioc(p', p", P) (A119) 
i?636(p',p",p) = -^{2FS36(p',p",p)-F56„(p',p",p) 

+F56,(p', p", p)} FB9c(p', P", P) (A120) 



5443(P',P",P) = ^FS6;.(p',p",p)m(p',p",p) (A121) 

5453(p', P", P) = -5354(p', P", P) = 2{(p' + p") . p}FB,{p', p", p) (A122) 

5463(p', P", P) = 5364(p', P", P) = -2{(p' - p") • p}FBr(p', p", p) (A123) 

i?553(p', P", P) = {2FS3„(p', p", p) - Fi?6c(p', P", P) 

+F56„(p', p", p)} FBship', p", p) (A124) 

i?563(p',p",p) = ^{2Fi33„(p',p",p) + Fi?6c(p',p",p) 

+FS6„(p', p", p)} FBiobip', P", P) (A125) 
5653(p', P", P) = ^ {2FS3„(p', p", p) - FBUp', P", P) 

-FBeaip', p", p)} FBnbip', p" , p) (A126) 

i?663(p',p",p) = -^{2Fi?3„(p',p",p)+Fi?6c(p',p",p) 

-FB,,{p', p", p)} FB,,{p', p", p) (A127) 

5334(p',p",p) = -^FS6„(p',p",p)FS6c(p',p",p) (A128) 

5335(P',P",P) = S353(P',P",P)=5633(P',P",P) 

= 5333(p', P", P) = 2FBr{p', p", p) (A129) 

5336(P',P",P) = B363(P',P",P)=S533(P',P",P) 

= -2FS7(p',p",p) (A130) 

i?343(p',p",p) = -^FS66(p,p",p)i^i?6a(p',p",p) (A131) 

5433(P',P",P) = -^FS6e(p',p",p)FS65(p',p",p) (A132) 
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5455(P', P", P) = {2F53c(p', p", P) - FBebip', p", p) 

-F56,(p',p",p)}' (A133) 

5456(P',P",P) = -^{2FS3c(p',p",p)+Fi?66(p',p",p) 

-FS6,(p',p",p)}' (A134) 

i?465(p', P", P) = {2F53,(p', p", P) - FB,,{p', p", p) 

+FS6e(p',p",p)}' (A135) 

5466(P',P",P) = -^{2FS3e(p',p",p)+FB66(p',p",p) 

+Fi?6c(p',p",p)}' (A136) 

i?545(p , P", P) = -4 {2^i?36(p', P", P) + FBUp', P", P) 

Id 

-F56,(p',p",p)}' (A137) 

5546(P',P",P) = -^{2FS3,(p',p",p) + FS6„(p',p",p) 

+FS66(p',p",p)}' (A138) 
5645(P', P", P) = {2FB,,ip', p", p) - FB,a{p', p", p) 

-F56,(p',p",p)}' (A139) 

5646(p', P", P) = {2FS3,(p', p", p) - FBUp', p", p) 
Id 

+FS6;,(P',P",P)}' (A140) 

i?554(p',p",p) = -^{2F53a(p',p",p)-Fi?6c(p',p",p) 

+FS6„(p',p",p)}' (A141) 

5564(P',P",P) = --^{2FS3„(p',p",p)+F56e(p',p",p) 

Id 

+FBUp',p",p)}' (A142) 

i?654(p',p",p) = -^{2F53a(p',p",p)-Fi?6c(p',p",p) 

-FS6„(p',p",p)}' (A143) 

5664(P',P",P) = -^{2Fi?3a(p',p",p)+i^S6c(p',p",p) 

-F56„(p',p",p)}' (A144) 

i3445(p', p", P) = -{(p" + P) • p'VFBrip', p", p) (A145) 

i?446(p',p",p) = -{(p" - P) ■ p'}^m(p',p",p) (A146) 

5454(p', p", P) = -{(p' + P") • pj'FB^ip', p", p) (A147) 

5464(p', P", P) = -{(p' - P") • pVFB^ip', p", p) (A148) 
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-B544(P',P",P) 
^644(P',P",P) 

5444(p',p",p) 
^566(P',P",P) 



-{(p +p)-p"rm(p',p",p) 

-{(p'-p)-p'rm(p',p",p) 

^656 (P', P", P) = 5665(P', P", P) 
^555(P',P",P) = -4-F^7(p',P",P) 



(A149) 

(A150) 
(A151) 

(A152) 



Appendix B: Coefficients for the deuteron 

In this appendix we present the expressions A'^ and B'^ given in Eqs. fl2.12p for the 
deuteron. 

Kip) = 3 
Kip) = <(p) = o 

AiM = (B1) 

The coefficients Bf--,{p,p') are exphcitly calculated as 



5m(P,P') 


= ^121 ( 


p,p 




= 




bU{p,p') 


= (p X 




BLiP,P) 


= (p'+ 


p)^ 


bUp^p') 


= (p'- 


p)^ 



BUp, p') = BUP; P ) = 4(p ■ p')' - Ip'p" 

i?232(P,P') = -8P-P (pxp')' 

BUp,p') = -ypV-4(p-pr + fpV(p-pr 

bUp, pO = I pW + P") - fp'p'\p ■ p') -lip'+ p")ip ■ P'f 

BUp. p) = \ pW + P') + ypV(p ■ P ) -\{P' + P"){P ■ Pf (B3) 



i?2n(P,P') = 52'2l(P,P') = i?23l(P,P0 = O 

i?2'4l(P,P') = -\P' {p'^Pf 

BUp, p) = Ip' + yP' (P • P') + 4(P • P')' - \p'p" 

BUp. p) = Ip' - fp' (P ■ P) + 4(p ■ pO' - Ip'p" (B4) 

The expressions for the functions -Bffc2) ^ — 6, can be obtained from the functions B2]^^ by 
replacing p ^ p'. 
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Appendix C: Example of a chiral potential 

For this particular example we will use the next-to-next-to-leading order (NNLO) chiral 
potential from Ref. {l2t]. 

The leading-order (LO) NN potential in the two-nucleon center-of-mass system (CMS) 
reads (Tsj: 

T. 1 ^i'^i-qcr2-q Cs ^ Ct 

- 'J^'Wi q2 + M2 • + (2^0^ + (2^0^^^^ ■ ' ^^^^ 

where q = p' — p, and m^, and denote the pion mass, the pion decay constant and 
the nucleon axial coupling constants. At next-to-leading order (NLO) a renormalization 
of the low energy constants (LECs) is required and the contribution from the Goldberger- 
Treiman discrepancy leads to a modified value of qa- The remaining contributions to the 
NN potential at this order are 



1 Ti-T2 



2/r- 4 A 2 i\ I 2/00 4 in 2 i\ , ^'^9a^ 



Vnlo = -J^,^^. L\q) \Ami{hg\ - Ag\ - 1) + c^\2^g\ - 10g\ - 1) + 4^, + 

L (g) CTi • q 0-2 ■ q - CTi ■ cr2 q 



(27r)3 647r2F4 

, ^1 2 , C'2 2 , / ^3 2 , ^4 , 2\ 

where g = |q| and k = ^(p' + p). The loop function L^{q) is defined in the spectral 



function regularization (SFR) as [12 

L\q) = eCA - 2m^) ^ In ' ' Z?'^^ ' (C3) 



a; , A^a;^ -f- 02^2 _^ 2kqujs 

— In 

2g Aml{JS? + 



with the following abbreviations: uj = \J 4m'^ + q ^ and s = y A^ — 4m^. Here, A denotes 
the ultraviolet cutoff in the mass spectrum of the two-pion-exchange potential. 

The contributions at NNLO again lead to the renormalization and/or redefinition of the 
LECs Cs, Ct, Ci,...Ci. The only new momentum dependence is due to the following 
terms: 

Vnnlo = --^,j^,{2ml{2c,-c,)-W){2ml + ci')A\q) 

- ^j||jTi-T2(4m^ + gV^(g)(^i-q^2-q-qVi.a2), (C4) 

where Ci, C3, C4 are new ttN LECs and the loop function A^{q) is given by 

A^{q) = 9{A - 2m^) arctan ~ ^"^-) . (C5) 

2g g2 _|_ 2Am7r 



29 



The expressions of the potential given in Eqs. ( ICip . f lC2p . and f lC4p show that this potential 
can be readily expressed in the operators Wj, j = 1,6 of Eq. fl2.2p . 
The chiral potential we consider in this example is the sum of 



y = Vlo + Vatlo + yNNLO- 



(C6) 



It requires regularization when inserted into the Lippmann-Schwinger equation, which is 
achieved by introducing a regulated potential of the form 



Ke.(p ',p) = e-(^'^/^^)y(p',p)e-(^^/^^) 
with the cut-off parameter A. 



(C7) 



Appendix D: Scalar functions for the Bonn B potential 

For the convenience of the reader we give the expressions for the Bonn B potential from 
Ref. 15|] in a form which is more suited for our three-dimensional calculations. The ex- 
pressions for the exchange of pseudo-scalar (ps), scalar (s), and vector (v) mesons are given 
by 



^ps(p',p) 
K(P,P) 

v;(p',p) 



Hps 



m m 




ps 



(27r)34m2 y E' y E {p' - p)^ + WW 



9t 



m 




(27r)34m2 y E' \ E {p' - p)2 + ml WW 

1 F^[ip' - pf\ {glO,, + 2g,,W,t + flOu) 

{27i)Hm^\l E'\l E {p' -py + ml WW 



(Dl) 



where rria are the masses of the exchanged mesons, m the nucleon mass, E = m? + p^ 
and W = m + E. The crucial quantities are the operators Ops, Os, O^^, O^t and Ou, which 
are given in terms of the Dirac spinors as 



Ops = 4m'W'W u{p')Yu{p)u{-p')Yu{-p), 

Os = -AmVW u{p')u{p) u{-p')u{-p), 
O,, = Am^W'W u{p')-i^u{p)u{-p')-i^u{-p). 



(D2) 

(D3) 
(D4) 



O^t = mW'W \4:mu{p')'y^u{p)u{-p')'yf,u{-p) 



- w(p')7''^^(p)^^(-p') 



u{-p) 



- HP 



(E' - E){g'^ - + {Pi + p'lT 



u{p)u{-p)-f^u{-p) 



(D5) 



Ott = W'W \4m^u{p'h^u{p)u{-p'hM-p) 



— 2m u(p')7^u(p)m(— p') 



iE'-E)ig'^~^,j') + ip,+p',] 



u{-p) 
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— 2mu{j>' 

+ ^(P') 
X u{—p 



w(p)m(-p')7mm(-p) 



i^E' ~E){g'^~ri') + {Vi+v'ir 



u{-p] 



(D6) 



with (pi + pi'Y = {E + E',p + p') and {p2 + p'^Y = {E + E', -p - p'). 

These operators act in the spin spaces of nucleons 1 and 2: the bihnear forms built with 
m(p') . . . u{p) contain cri as acting in the spin space of nucleon 1 and the bihnear forms with 
u[—p') . . .u{—p) contain (T2 and act in the spin space of nucleon 2. The spinors u{q) are 
normalized according to the definitions given in Ref. 23|] and explicitly given as 



Each vertex is multiplied with a form factor 




(D7) 



^a[(p'-p) 



A?. - ml 



IP' - Pj 



2n 



(D8) 



where the values of n and the cutoff parameters are given in Table [TTl 

Note that for the three iso- vector mesons (tt, 6 and p) contributing to the Bonn B poten- 
tial, expressions (D2)-(D6) are additionally multiplied by the isospin factor t(1) • t(2). 

In Ref. [il this potential was presented in a different operator form. However, in that work 
one of the six operators was chosen to be cri ■ (p + p') cr2 ■ (p' — p) + cri ■ (p' — p) cr2 ■ (p + p') , 
which is an operator that violates time reversal invariance. In practice, this operator is 
always multiplied with the term [p'"^ — p'^), which also violates time reversal invariance. 
Therefore, the entire term is invariant as it should be. In principle it is not desirable to work 
with symmetry violating operators, thus we prefer to use the operators from Eq. fl2.2l) and 
rewrite 

CTi ■ (P + P') cr2 ■ (p' - p) + o-i ■ (p' - p) 0-2 ■ (p + p') = 

-4(p X p'^2 ( 



p/2 _ p2 



(Tl (T2 + ,2 _ 2 ^1 ■ (P + P ) ^^2 ■ (P + P ) 



p' — p 

-(Ti ■ (P - P ) 0-2 ■ (p - p ) + ^2 2^1 ■ (P X P ) Cr2 ■ (P X p 



p' — p^ 



p — p^ 



(D9) 



which is an identity for (p + p')(p' — p) = p'^ — p"^ 7^ 0. Inserting Eq. (lD9p into the 
expressions given in cancels the factor p'^ — and one obtains the following expressions 
for the operators Oa from Eqs. flD2p - flD6l) in terms of the operators wj = Wj{(Ti, a2,p' ,p) 
from Eq. (Q: 



O 



ps 



1 + 



2m 



E' + E 



II \2 /2 2 
(P ■ P) -P P 



^2+1 + ^ 



2m 



E' + E 



1 I , ,9 / 2m 



[p"+p'-2{p'-p)] }W5 
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+ Wr + (^1 + \p" +p' + 2(p' . p)] (DIO) 



Os = -[WW - (p' • p)]^Wi - [WW - (p' • P)]W3 + W4 

Oyy = {[WW + (p' • p)]2 + w'Y + wy^ + 2iy'iy(p' ■p)]wi 



1 / 2m 
+ 2^ + ^ 



^ [p" + - 2{p' . pf]^w. 



-[^WW+{p'-p)]w,-{2 + ^^ 

-\ [-{W - Wf + (l + - 2(P' • P)]} 

{-(W^' + wf + (l + + 2(p' . p)]} 



(Dll) 



(D12) 



O 



vt 



wy + wy^ - ^' ^ (wy - wy) 

+ TW'W [2WW + (p' • p)] - ^' \W^W'' - (p' • p)'l \ wi 

m ' 

+ |2iy'W^(p' • p) - ^ {W^ + W^) {p'^+p^) 

W + W 

2W'W + (p^ ■ p) 



+ 



[W^ + W^- 2m{W + W)] \ Wi 



2m{E' + E) 

-\ [-{W - wf + (l + b" + - 2(P' • P)] 

- [ly' V + y - (w^" + vp^') (p' ■ p)] j 



\ i^-{W + wf + [p'^+p' + 2(p' . p)] 



1 



m(E' + E) 



[wy + wy^ + + w^) (p' • p)] \ w. 



(D13) 



[W^'H^ + (p' • p)]' + 2(2- ^'^^ ^ [W^W^ - (p' • pf] 
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+ <^3 + 



'i[W'W - m{W' + W)] + (p' • p) 



2m^ 



+ 



1 + 



{W - W) 



21 



Am? 



{W'Y + W^p'"^) + 2 



[ly'w-(p'-p)] 
^ _ {W - wf 

Aw? 



W'WijI ■ p) 



w -w 



m 



+< 



_ {w - wy 

Am? 

^ I {W - W) 
Am^ 



W'W{^' ■ p) - 
m2 + E'E 



1 + 



{W - Wf 



E' + E 



Am^ 

+ m ) {W'p^ + Wp'"^) 



2m 



E' + E 



IP ■ P) 



m{E' + E) 



[W'Y + wY - {w' + w^) (p' • p)2] 



'W2 



+{ -iy'iy-(p'-p) + 2 2- 



+ <^3 + 



W' + W 



m 



(P' ■ P) - 2 



^ _ jW' - Wf 
Am^ 



W'W 



3[iyW - m{W' + W)] + (p' • p) 
2m2 



[ly'iy - (p' • p)] !>«;3 

W' + w' 



2m'^ 



m 



+ 



1 + 



+ 



+ - 



2 

1 
2 

1 
2 

1 



1 + 



(W^' - Wf 
Am^ 

_ iW' - Wf 
Am^ 

iw' -wy 

Am^ 

iW' -Wf^ 



2m 



1 



W'W 



E' + E) m{E' + E) 
1 



m{E' + E) 
2m 



[wy + wY - {w'^ + w^) (p' ■ p)] 



1 - 



1 + 



4m^ 

{W -W)^' 
Am^ 



E' + E 
W'W - 



+ p'- 2(p' ■ p)] - W'^ - W^ 

1 



m{E' + E) 



[W^p^ + w^p'^ + {w'^ + w^){p' ■ p 



1 + 



2m 



E' + E 



[p" + p' + 2(p' ■ p)] - W" - W- 



-Wq. (D14) 



The values of the parameters are given in Table [TTl 
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